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The breaking of Lorentz symmetry via a dynamical mechanism, with a tensor field which takes 
on a non-zero expectation value in vacuum, has been a subject of significant research activity in 
recent years, fn certain models of this type, the perturbations of the "Lorentz-violating field" 
about this background may be identified with known forces. I present the results of applying this 
interpretation to the "generalized bumblebee models" found in a prior work. In this model, the 
perturbations of a Lorentz-violating vector field can be interpreted as a photon field. However, the 
speed of propagation of this "bumblebee photon" is direction-dependent and differs from the limiting 
speed of conventional matter, leading to measurable physical effects. Bounds on the parameters of 
this theory can then be derived from resonator experiments, accelerator physics, and cosmic ray 
observations. 
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I. INTRODUCTION 

The experimental signatures of a violation of Lorentz 
symmetry have been extensively sought for in recent 
years (see [T] and references therein.) The primary 
paradigm for examining the physics of such effects is 
the "Standard Model Extension" (SME) HIS]. Broadly 
speaking, in the usual picture of the Standard Model, 
one writes down a list of field combinations that are 
renormalizable and invariant under Lorentz symme- 
try (as well as under various other desired symmetries 
amongst the fields), assigns a coefficient to each one, and 
writes down the Lagrangian as the most general linear 
combination of these terms. The values of these coeffi- 
cients are then to be established by experimental mea- 
surements. The SME "extends" this paradigm by relax- 
ing the requirement that the field combinations in the 
Lagrangian be Lorentz- invariant. Since the Lagrangian 
itself should still be a Lorentz scalar, the coefficients 
of these new terms must have non-trivial tensor struc- 
ture (rather than being Lorentz scalars as in the original 
Standard Model.) These new coefficients (or, more ac- 
curately, their components in some reference frame) can 
then in principle be measured via experiment. 

While this method works well for the purposes of par- 
ticle physics, it becomes somewhat problematic when 
we attempt to extend it to gravity. With a flat metric, 
it is legitimate to view the new Lorentz-tensor coeffi- 
cients in the SME as constants throughout spacetime. 
The notion of a constant tensor field on a flat back- 
ground spacetime is well-defined; we simply require that 
(for example) daV^ — 0. However, once we allow for a 
curved background, it is no longer so simple to find a co- 
variantly constant non-zero vector field (i.e.. Vat''' = 0); 
indeed, such a vector field may not even exist on an arbi- 
trary curved background. Moreover, general arguments 
involving the Bianchi identities [1] imply that any "back- 
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ground tensor field" that couples directly to the curva- 
ture in the Lagrangian must satisfy certain differential 
conditions; we cannot simply write down some fixed ten- 
sor fields on our manifold and proceed from there. 

The standard way to solve the problems arising in 
making the metric dynamical is to also promote the SME 
coefficients to dynamical fields. As the Lorentz-violating 
fields are now dynamical, there is no reason to expect 
them to be covariantly constant, and the geometric con- 
sistency of the field configuration with the Bianchi iden- 
tities is automatic. One then constructs the theory such 
that these fields take on some non-zero value in the limit 
of no conventional matter and flat spacetime (hence "vi- 
olating" Lorentz symmetry by taking on a non-invariant 
background value.) The flat-spacetime SME is recov- 
ered by constructing an effective field theory about this 
background, where the metric and the Lorentz-violating 
fields are held fixed but the other fields in the theory are 
allowed to vary. 

While this promotion of Lorentz-violating coefficients 
to Lorentz-violating fields solves the above problems, it 
does require some care. In particular, the requirement 
that the dynamics of a Lorentz-violating vector field not 
fundamentally change the dynamics of the metric re- 
stricts us to a small subclass of all conceivable vector 
models [5] . The resulting models have the property that 
the dynamics of the field are decoupled from the dynam- 
ics of the metric (at least at the linearized level.) 

This decoupling between the Lorentz-violating field 
and the metric had been previously seen in simpler vec- 
tor models known as "bumblebee models" |S1 [7] • It was 
further noted that the linearized equations of such a sys- 
tem were (under minor auxiliary conditions) precisely 
those of linearized Einstein-Maxwell theory. The pertur- 
bations of the Lorentz-violating vector field could then 
be interpreted as the photon field in such a theory. Un- 
der such an interpretation, the dynamics of the known 
long-range forces would be the same as in conventional 
Einstein-Maxwell theory, up to small Planck-suppressed 
deviations. 

The class of models found in '5] included the 
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previously-known bumblebee models, and so were 
dubbed "generalized bumblebee models". However, it 
was noted in that paper that attempting to extend the 
"bumblebee photon" interpretation to these generalized 
models could lead to readily observable effects. Specifi- 
cally, in a generalized bumblebee model the metric per- 
turbations and conventional matter will "see" a different 
metric than the photon field will; in other words, the 
"speed of light" would differ from the "speed of gravity" 
and the limiting speed of conventional matter. 

The present work elaborates on the above speculation. 
Specifically, we will derive the observable consequences 
of the "generalized bumblebee photon" theory, and place 
bounds on the parameters of the underlying Lorentz- 
violating vector field. Section [Tl] reviews the derivation 
of generalized bumblebee models and describes the "pho- 



ton interpretation" mentioned above. Section III derives 
the possible experimental signatures for such theories, 
both in the context of the SME and in terms of par- 
ticle kinematics. Finally, Section |IV| examines the cur- 
rent experimental and observational bounds on Lorentz- 
symmetry violation in the photon sector, and derives 
bounds on the parameters of the underlying Lorentz- 
violating field. 

We will use units in which h — 8ttG = I through- 
out. We will also set c = 1 in Section |llj however, to 
avoid confusion between the various limiting speeds in 
subsequent sections, we will explicitly include all such 
speeds in our equations. Sign conventions concerning 
the metric and the curvature tensors are those of Wald 
[5], except in the Appendix where we use the signature 
(+,-,-,-)• 



II. GENERALIZED BUMBLEBEE MODELS 

If we limit ourselves to theories of second differential 
order containing a single vector field Ba along with the 
metric gab, the most general model of dynamical Lorentz 
symmetry breaking has the action 

I d^X^iR + J'^'^'^Rabcd 

+ IC''''''VaBbVcBd-V{B^)) (I) 

where R is the Ricci scalar derived from gab, Rabcd the 
Riemann tensor, and J"''^"^ and JC"'^'^'^ are arbitrary ten- 
sors constructed locally out of Ba and gab- The potential 
V^B"^) is taken to vanish and to be minimized at some 
non-zero value of its argument. Under these assump- 
tions, any field configuration with 



gab = riab and Ba^ Ba, 



(2) 



where Ba is a constant non-zero vector field with 
V{B^) = Vi^B"^) = 0, is the "natural" solution of the 
equations of motion. This non-zero vector field then 
provides a "preferred direction" in spacetime. 



While a Lagrangian of the form ([T]) is indeed the most 
general form for the Lagrangian, it was shown [5] the 
equations of motion derived from a completely arbitrary 
Lagrangian have certain less-than-dcsirable properties. 
In particular, when varying the kinetic term for the vec- 
tor field K.'^^'^'^V aBbV cBd, we find that it gives rise to 
second derivatives of Ba in the Einstein equation, and 
that these terms cannot in general be eliminated via the 
vector equation of motion. We are thus left with a sit- 
uation in which the dynamics of the vector field are in- 
herently coupled to those of the metric. 

However, a certain class of vector models do not ex- 
hibit this coupling. In particular, if the kinetic term for 
the vector field is of the form 



IC'^'^'^WaBbWcBd = 



'Cg'^^r^FacFb^ 



(3) 



where 



~ab 



(4) 



and Fab = 2V[ai?f,], then the equations of motion for 
the vector field and the metric decouple.^ (The nega- 
tive sign is chosen for agreement with convention; C may 
be positive or negative.) This decoupling is not terribly 
surprising when one remembers that the field strength 
Fab is proportional to the exterior derivative of a one- 
form, and thus is independent of the derivative operator; 
varying the metric (and its associated covariant deriva- 
tive operator) therefore does not give rise to any terms 
containing the second derivatives of the vector field as 
it does in the general case. (See §4 of [9] for further 
discussion.) 

For the remainder of the paper, we will restrict our at- 
tention to theories with "pseudo-Maxwell" kinetic terms 
of this type. We will also take J°-^^'^ to vanish; the 
primary effect of such terms in pseudo-Maxwell vector 
theories is to modify the "effective Einstein equation" 
[SI [TO] , but gravitational effects are not the primary fo- 
cus of this paper. Our models will thus be those derived 
from an action of the form 



-g {R - Cr'r^FacFbd - V{B^)) . (5) 



^ Although we will take C, and /3 to be constants for most of the 
paper, it is possible that they might themselves be functions of 
. If this is the case, the decoupling still holds at the level of 
the linearized equations, with C, and /3 being replaced by C,{B^) 
and P{B'^). 
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A. Linearized equations 

Varying g"'' and Ba in the action ^ , we find that the 
full equations of motion are of the form 



Gab = 2C( FacFM^'' " -gabFcdFefg^'g''^ 



2PB^aFb)cFderB^ 



and 



V'iB^)BaBh- -gabViB^) (6) 



PF\F,d~g'"'B'' = ^^V'{B^)B\ (7) 



We now linearize these equations about the back- 
ground described above, i.e., we let 



and 



gab — Vab + hab 



Ba — Ba + Aa, 



(8) 



(9) 



where Ba is a constant vector field on Minkowski space- 
time and hab and Aa are considered to be "small" . Re- 
quiring that this field configuration be a solution when 
hab and Aa vanish implies that V{B'^) = V'{B^) = 0, as 
noted above. If we then linearize the equation Q about 
this background, we obtain the linearized Einstein equa- 
tion 



SGab = V"iB^)BaB,6iB^), 



(10) 



where SGab is the linearized Einstein tensor (in terms of 
derivatives of hab) and S{B^) is the linearized variation 
in the norm of B„ , 



6{B^) = - B''B''hab 



Linearizing the vector equation of motion ([t]), mean- 
while, yields 

r^'rf^dbid.Ad - ddA,) = ^V"{B^)B-S{B^), (12) 

where the background value of g"'' Q , the "effec- 

tive metric" for 



~ab ab 



a r>b 



PB^B 



(13) 



B. Charged-dust equivalence 



The linearized equations of motion (10 1 and ( |12| , 
though simpler than the full equations (|6| and ([t]), are 
still somewhat complex. To gain some intuition about 



their solutions, let us first consider the case of the the- 
ory in which /? = 0; in this case, the linearized vector 
equation becomes 



d'ha 



1 

2C 



V"{B^)Ba5{B^) 



(14) 



where fab = 2i9[o^h]. It was noted by Bluhm, Fung, 
and Kostelecky [SI [7] that there is a one-to-one cor- 
respondence between solutions of these equations with 
S{B^) = (in a particular gauge) and solutions of con- 
ventional linearized Einstein-Maxwell theory (in a par- 
ticular gauge.) Specifically, if we apply an infinitesi- 
mal diffeomorphism (parametrized by a vector field ^°') 
to our background field configuration gab = Vab and 
Ba = Ba, these fields transform as 



Vab ^ Vab 



and 



Ba ^ Ba + B^da^b, 

since 9a ^6 = 0. The transformations 

hab — >■ ^Ifc = hab + '^d^aCb) 
Aa^A'^^Aa+ B'daCb 



(15a) 

(15b) 

(16a) 
(16b) 



are therefore gauge transformations and do not affect 
any physical quantities. In particular, for a solution of 
(10 1 and (14) with S{B^) — 0, we can apply a gauge 
transformation with chosen such that 



B^di^a^t) = -^B'^hab, 



(17) 



thereby putting hab in the axial gauge (i.e., B'^h'^i^ = 0.) 
Importantly, under such a gauge transformation, we will 
also have 



(11) B-A'^ = B-Aa 



B-B^da^b = B'^Aa - ^B-B'hab = 0, 

(18) 



since we are assuming that 6{B ), as given in ( [Tl] ), 
vanishes. Thus, putting hab in axial gauge automati- 
cally also puts Aa in axial gauge if 6{B^) = 0. More- 



over, in this case the equations of motion ( 10 ) and ( 14 ) 



are simply the source-free Einstein and Maxwell equa- 
tions. Thus, every solution of ( 10 1 and ( 14 ) for which 



S{B ) = can be mapped to a solution of conventional 
Einstein-Maxwell theory for which both the metric per- 
turbation and the vector field are in axial gauge. This 
mapping can also be seen to go the other way (again up 
to gauge transformations): given a solution of conven- 
tional source-free Einstein-Maxwell theory, apply gauge 
transformations to both hab and Aa such that they are 
in axial gauge with respect to Ba- This gauge trans- 
formation guarantees that S{B^) vanishes, and thus this 
field configuration is also a solution of ( 10 ) and ( 14 ) with 
(5(^2) =0. 
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As it turns out, this mapping can be extended to the 
case where S{B'^) ^ 0, at least in the case where Ba is 
timehke. Let us suggestively define^ 



2< 



and 



Ua=± 



Ba 



(19) 



(20) 



(21) 



The signs of pe and Ua are chosen to be positive if i?° 
is future-directed and negative if it is past-directed; in 
other words, u° is defined to be future-directed. Rewrit- 



tions become 



ing ( 10 1 and ( 14 ) in terms of these quantities, the equa- 



SGab = PmUaUb 



PeUa 



(22a) 



(22b) 



which are easily recognizable as the equations of motion 
for the perturbed metric and the vector field Aa in the 
presence of charged dust. By construction, Ua is a unit, 
future-directed timelike vector. The charge-to-mass ra- 
tio of the dust is constant, and is given by 



Pe/Pn 



1 



2C\/=i?2' 



(23) 



By applying the Bianchi identity to (14 1, we obtain 

B'^daiSiB')) - 0, (24) 

which guarantees that pm and pe are constants along the 
worldlines parametrized by u"". 

We therefore conclude that any fields hat and Aa sat- 
isfying ( 10 1 and ( 14 1 can be mapped to a solution of 



conventional Einstein-Maxwell theory with a charged 
dust source, where the dust moves along the worldlines 
parametrized by B"^, and its mass density and charge 
density are given by (19) and (20 1 respectively. As in 



the case of vanishing T(i3^), this correspondence goes 
the other way as well. Suppose we have a solution 
{hab,Aa} of the linearized Einstein-Maxwell equations 
with a charged-dust source with mass density and 
charge density pe, with p^ and pe satisfying (23). We 



can perform a gauge transformation on the Maxwell 
field, Aa ^ Aa + 9aA, with A satisfying 



B^'daX 



Pn 



2V"{B^)B^ 



B'^Aa 



1 



B'^B'^ha 



(25) 



(This does not uniquely determine A, of course, but we 
only require A to exist.) Under this gauge transforma- 
tion, the fields hab and Aa will satisfy 



S{B' 



2B''A„ 



B^'B'^hn 



V"{B^)B^' 



(26) 



We can then see that in this gauge the fields hab and 
Aa satisfy our original equations (10) and (14). 



This 



le original 



correspondence is easily seen to agree with t : 
correspondence [Hl[7j in the case where S{B^) = 0. 

This correspondence, between solutions of our lin- 
earized equations ([lO|) and ( 14 1 and those of Einstein- 



Maxwell-charged-dust systems, can be then used to gain 
some intuition about the behaviour of our system."^ In 
particular, this correspondence justifies the tactic (used 
in [HIT]) of simply setting the (5(S^) term to zero. One 
might have been concerned that this set of solutions 
was unstable, in the sense that a solution with 5{B^) 
initially small but non-zero might evolve to a solution 
with large S{B^). This new correspondence shows that 
this is not the case, since sufficiently small S{B^) on 
the "bumblebee" side corresponds to small sources on 
the Einstein-Maxwell side, and a solution of the conven- 
tional Einstein-Maxwell equations with a small source 
will be "close" (in an appropriate sense) to a solution of 
the Einstein-Maxwell equations with no sources. From 
here on, we will assume that 6{B^) is negligible unless 
otherwise stated. 

A similar correspondence was noted by Jacobson and 
Mattingly [TT] in their studies of "Einstein-aether the- 
ory." In this case, however, the vector field they were 
examining served a dual purpose as both the vector po- 
tential and the dust worldlines; this implied, in par- 
ticular, that their "dust" was dynamical rather than a 
fixed background source. Since these two vectors are 
not in general aligned in an arbitrary Einstein-Maxwell- 
charged dust system, the correspondence found in [11] 
was therefore not one-to-one (even after taking gauge 
transformations into account.) Since our correspondence 
uses the fixed background as a "source" for the lin- 
earized perturbations, it does not run into this difficulty; 
any solution of the linearized bumblebee equations can 
be gauge-transformed into a solution of the Einstein- 
Maxwell equations with a charged-dust source, and vice 
versa. 

Finally, recall that all of the analysis in this subsection 
has been done assuming that the constant /3 vanishes. 
The above analysis changes in two main ways if /3 ^ 
0, one less important and one more important. The 
first is that the charge-to-mass ratio of the dust in the 
above correspondence changes. The linearized Maxwell 



^ Note that the quantities pm and pe are proportional to the quan- 
tity 13 defined in Eqn. (68) of [7]. 



Our correspondence also seems to work in the case of spacefike 
Ba- However, in this case the "dust" sources will be moving 
along spacelike worldlines, a situation of which it is less common 
to have an intuitional understanding. 
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equation ( 12 ) in this case becomes 



2C(1 + /3S2) 



(27) 



with Ua defined as in (21 1. (To see this, multiply by 
the tensor fjab defined such that fjabrf"^ — Sa'^.) Thus, 
the "charge density" defined in ([20| is multiplied by a 
factor of (l+/3-B^)^^ when we pass to the general case of 
non- vanishing (3.^ The rest of the above argument holds, 
however; in particular, we are still justified in assuming 
(5(5^) to be negligible. 

More importantly, however, when /3 7^ the vector 
perturbations Aa will not propagate with the same ve- 
locity as those of the metric. Instead, the metric pertur- 
bations will propagate along the light-cones of the usual 
flat metric 77°'', while the vector perturbations will prop- 
agate along the light-cones of the "bumblebee metric" 
77°^. Assuming that any matter sources are minimally 
coupled to the "Einstein metric" gab used in the action 
([5]), and that their kinetic terms are not directly cou- 
pled to B°-, this also imphes that the limiting speed of 
conventional matter will be different from the limiting 
speed of the vector perturbations. The observational 
consequences of this fact will be explored in the next 
section. 



III. LORENTZ- VIOLATING PHOTONS 



A. Bumblebee photon theories 



effects of the Lorentz-violating field might still be 
observable via gravitational effects |10] , but would 
not give rise to forces between particles of "con- 
ventional" matter. 

• The Nambu-Goldstone modes are "eaten" by an- 
other field via a Higgs mechanism. This turns out 
to be impossible [12] in the context of spacetime 
with a Riemann metric, though it is possible in 
Riemann-Cartan spacetimes with a dynamical tor- 
sion field [6j. We will not consider this possibility 
further here. 

• The massless field Aa couples directly to conven- 
tional matter, giving rise to a long-range "fifth 
force". For example, the field A'^ could conceivably 
couple to leptons but not quarks (or vice versa), 
it could couple differently to first-generation par- 
ticles than to second-generation particles; it could 
couple only to strange quarks; and so on. The 
large number of experimental signatures that con- 
ceivably could arise in such scenarios are, unfor- 
tunately, outside the scope of this paper; models 
along these lines have been explored in [T31 [T3] . 

• The massless field Aa couples directly to conven- 
tional matter, and can be identified with a known 
force. The obvious candidate here (as may have 
been telegraphed by the choice of notation) would 
be the photon field [51 \7\ . This interpretation will 
be the focus of the rest of this work. 



We found in the last section that linearized solutions 
of the equations (|6| and ([T]) (about a background where 
the metric is fiat and the vector is non-zero) can be taken 
to satisfy the equations 



SGab - 



and 



v'^dbfc, 



0. 



(28a) 



(28b) 



From the perspective of particle physics, these are mass- 
less fields (or more precisely, Nambu-Goldstone modes 
arising from a spontaneously broken symmetry.) One 
can envision a number of distinct possibilities concern- 
ing the effects of such fields on the theory: 

• The field Ba does not directly couple to conven- 
tional matter. In this case, we would not have 
detected its effects in particle experiments. The 



* We are of course assuming here that /3_B^ ^ —I. In the case 
where I3B^ = —1, the inverse metric defined in |13| becomes 
degenerate, and | |12[ l cannot be viewed as an evolution equation. 
We will assume hereafter that /3 and are chosen such that 
r)"'^ and j)"' have the same signature. 



One might ask whether it is self-consistent to demand 
that ry"'' serve as the "conventional matter metric" while 
simultaneously requiring that the bumblebee perturba- 
tions Aa serve as the photon. This self-consistency can 
be shown by examining the possible couplings between 
B'^ and the fermion fields in the theory. Suppose we 
have a fermion field tp appearing in the Lagrangian with 
its standard kinetic term and with Ba coupling to its 
current: 

= '^h^dai^ + qBai^r^ 

= '-^p-f^'daij + qBai^l^lp + qAai^^lp. (29) 

We can then see that the decomposition of Ba into a 
background field plus a perturbation (identified as the 
photon) leads to the usual interaction term qAa'>jjj°"ip be- 
tween the bumblebee photon and the fermion. The sec- 



ond term on the right-hand side of (29 1, meanwhile, can 



be interpreted in the language of the Standard Model 
Extension (SME) 3J as a Lorentz-violating coefficient 
Oq = qBa. Through a redefinition of the spinor phases, 
the coefficients Oa in the SME can be made to vanish in 



fiat spacetime [2]. Thus, a term of the form (29) would 
give rise to a conventional photon-fermion interaction, 
without other observable effects in the fermion sector. 
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We could also envision having the fermion interact 
with the bumblebee field via a derivative interaction: 

= l^l'da^ + '-qcB''B''{l,^adb^J (30) 

where Qc is a coupling coefficient. In the language of the 
SME, such a term would give rise to a Cab coefficient for 
the fermion field ijj-^ It is precisely such a Cab term that 
would cause the "effective fermion metric" to differ from 
j^ab G j£ consider the bumblebee field as taking on its 
fixed background value, the above Lagrangian (30) can 
be rewritten as 



where fj'^'' = 77"'' + qcB°-B^. We could equally well de- 
fine fj°'^ to be our "fundamental metric" instead of 77"''; 
this essentially amounts to a rescaling of the coordinates 
[T5] . We would thus have a theory in which the electrons 
propagate with respect to the 
the photons propagate with respect to 



"fundamental metric" r;"'', 



~ab ^ ^ab _^ _ q,)B''B'' 



(32) 



and the metric perturbations propagate with respect to 



ah 



ab 



QcB'^B' 



(33) 



We can then see that a theory with a non- vanishing is 
physically equivalent to a theory with gc — J' 0, /3 — > /3 — 
Qc, and a "distorted metric" for the metric perturbations. 
As the remainder of the paper will not be concerned with 
the metric perturbations, we will therefore assume that 
Qc has been set to zero in this way, and we will use ry"^ 
to denote the "matter metric" . 



B. SME coefficients 

If Aa is to be interpreted as the photon field in our 
theory, we immediately note an important experimental 
consequence of this fact: the photon does not propagate 
along the null cones of the conventional matter metric 
ry"^, but rather along those of the distorted metric 77"^ 
defined in (13). This distortion will, in principle, be ex- 



perimentally detectable. The potential effects of a back- 
ground geometric structure on the propagation of pho- 
tons were explored in detail by Kostelecky and Mewes 
|16j . One starts with a photon Lagrangian of the form 

Id (34) 



C = FabF'^" + {k^r'^^FabF, 



^ In principle, we could also couple Ba to the axial fermion current 

_ 

'0757'''!/i or to a term of the form i'lp-fsfa db''P\ such terms would 
give rise to ba and dab coefficients in the SME, respectively. In 
this work, we will assume these vanish. 
® Such a term would also give rise to momentum-dependent 
fermion-fermion-photon vertices, as well as two-photon-two- 
fermion vertices; however, such terms would be nonrenormaliz- 
able, and therefore would be highly suppressed at low energies. 



(up to an overall normalization), where {kp)"'^^'^ ~ 
(/cp)!"''!!'^''! is symmetric under the exchange {ab} ^ 
{cd} and has vanishing double trace (i.e., {kp)°-''ab = 0.) 
The tensor (kp)'^'"^''' can then be decomposed into vari- 
ous "electric" and "magnetic" parts that determine the 
electric and magnetic susceptibility of free space, as well 
as vacuum birefringence effects. In our case, the effective 
fiat-space Lagrangian for Aa is given by 



C = r^v'^FabF^d 
which corresponds to a [kp)"'^'^'^ tensor of 



(35) 



(31) (kp) 



abed 



1 



2 ^ 



B^ 



(36) 

(The factor in the denominator arises from factoring out 
the overall normalization mentioned above.) For the 
purposes of comparison with experiment, the compo- 
nents of {kp)'^'"^'^ can be decomposed into four spatial 
matrices ke± and ko± and a trace component Ktr, as de- 
fined in Section II B of Ti)].^ In the current case, these 
work out to be 



0, 



(ke 



and 



V'J - p [ B'B^ - -5'^B' 



2 v(^°)' + r' 



(37a) 



(37b) 



(37c) 



(37d) 



where B denotes the spatial components of S", and we 
have defined 



1 + lB^- 



(38) 



We can then use experimental measurements (see [T] and 
references therein) of the components of ke± , ko± , and 
ktr to place bounds on the parameters /3 and B'^ of our 
theory. 



C. Particle kinematics 

In a theory in which the limiting speed of a charged 
particle species is identical to the speed of the photon, 



Note that we have also implicitly defined a reference frame in 
defining these as "spatial" matrices. In what follows, we take 
this frame to be the standard Sun-centred frame, where the Sun 
is at rest, the Z-axis points towards the North Celestial Pole, 
and the X-axis points towards the Vernal Equinox. 
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it is kinematically forbidden for a photon to decay to 
that particle and its antiparticle, or for the charged par- 
ticle to radiate a photon. When the limiting speed of a 
charged particle differs from the speed of light propaga- 
tion, however, such processes are kinematically allowed 
(see Figure [ij) More precisely, if the speed of light 
in a given direction is greater than the limiting particle 
speed Cp in that direction, then photons above a cer- 
tain energy can decay; if in a given direction is lesser 
than Cp in that direction, then particles above a certain 
energy will undergo vacuum Cerenkov radiation. 

It is important to note that each of these processes 
is sensitive to /3 values of only one sign. The directions 
ka of bumblebee photon propagation are those for which 
fi"'''kaki, — 0. Denoting fc^ = {uj/cp,k), where Cp is the 
limiting speed of the particle species, we see that this 
"null condition" is equivalent to 



C.2=c2(fc2+/J(B-fc,)2) 



(39) 



This implies that if > 0, we will have = a;/|fc| > Cp, 
and photons of sufficiently high energy can decay to 
charged particles and anti-particles. Similarly, if /? < 0, 
we will have Cj < Cp, and charged particles of sufficiently 
high energy will lose energy to vacuum Cerenkov radia- 
tion.® 

What are the threshold energies for these processes? 
Denoting k = kk, where fc is a unit vector, we find that 



the photon dispersion relation ( 39 1 above can written as 



— Cpk 



I3B"{B ■ k) 



l + /3(-(S0)2 
1-/3(^0)2 



{B-kf 



(40) 

In the case of photon decay, the threshold photon energy 
will be that for which a photon with four-momentum 
hcpka can decay into a particle-antiparticle pair, with 
each particle having four-momentum pa = -^hcpka and 
rest mass m. Taking the norm of Pa with respect to 77"'', 
we find that the four-momentum of the each particle 
must satisfy 

PaPbr' =P''Pa+(3{B''Pa? - -m^ ^ + /3 {B^af . (41) 



The left-hand side of this equation vanishes (since 
PaPbff'^ oc kak}jf]°'^). Rewriting the right-hand side in 
terms of fc^, and defining rt" = yJ\P\B°- , we find that at 



photon 
mass shell ,. 



particle 
mass shell 




(a) 



particle 
mass shell 



photon 
< mass shell 




(b) 



FIG. 1: Mass shells and momentum vectors for one-to-two 
processes when the limiting particle speed Cp differs from the 
speed of light c-y = E^/p-^. (a) When > Cp (or, equiva- 
lently, /3 > 0), a photon (blue) can decay into two massive 
charged particles (red), (b) When c-y < Cp (/3 < 0), a charged 
particle with sufSciently high energy (red) can decay to a 
charged particle of lower energy (red) and a photon (blue). 



threshold 



The generahzed bumblebee photon model has what might be 
called a "homogeneous" photon dispersion relation: for all 
A > 0, if A;" is a valid four-momentum for a propagating photon, 
then so is Afc". Our discussion below can easily be extended to 
any model in which this is the case. For models in which this 
does not hold (see, for example, [17) V the "geometric" argu- 
ments used below to find the vacuum Cerenkov threshold can 
be adapted to analyze both photon decay and vacuum Cerenkov 
processes. 



{n^uj + [n ■ k)cpkY = 



2mCp 



(42) 



Applying the dispersion relation ( 40 1 then yields a 
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threshold energy of 



n° - (n- k)\Jl - {n°Y + (n • 



(43) 

In the hmit of the components of 71^ being much less 
than unity, this simplifies to 



2mci 



\TnP + fi ■ k\ 



(44) 



In the case of vacuum Cerenkov radiation, the thresh- 
old energy is best found using a geometric argument 
(see Figure [1]) Since momentum is conserved, the four- 
momentum of the photon will "connect" two points on 
the charged-particle mass shell. The threshold energy 
EcQ for vacuum Cerenkov radiation is therefore that 
point on the mass shell at which the slope of the tangent 
line {dEp/dpp) equals the slope of the photon mass shell: 
at any energy on the charged-particle mass shell with 
Ei > Eco, we can draw a secant line with the same slope 
as the photon mass shell that will intersect the charged- 
particle mass shell at a lower energy Ef < Ei. This will 
not, however, be possible for Ei < E^o- Performing this 
calculation, we find that the threshold energy Eco for a 
particle of mass m is given by 



Ecoik) 



1-C2(fc)/C2 



(45) 



where the (direction-dependent) speed of light c-y{k) is 
given dividing a; by fc in (40). Plugging this in, we find 



that the threshold energy is given by 



Ecoik) ^ mcl 



(n-k)- nO^l - (n0)2 + {ft ■ ky 



(nO)2 - (n ■ ky 



(46) 

where we have again used the rescaled vector = 
^y\P\B°-. If we again take the limit of small n^, this 
reduces to 



Ecoik) 



I'nP + n ■ k\ 



(47) 



In the case /3 < 0, the detection of a charged particle 
with energy Ep and mass m propagating in the direction 
k implies that the Cerenkov threshold energy Eco{k) for 
that direction is greater than Ep. Thus, we can say that 
such a detection restricts the components of n° to lie in 
the region 



+ n-k\ < 



E„ 



(48) 



This region is a thickened plane in n^-space, with a to- 
tal thickness of 2mc^/Ep. Multiple detections of charged 



particles coming from different directions k can then con- 
strain the parameters of our theory to a finite region of 
n^-space. Note that the quantity appearing on the right- 
hand side of ( [48| is simply the inverse of boost factor 7 
of the charged particle detected. 

Similarly, when /3 > 0, the detection of a photon with 
energy E~f propagating in the direction k implies that 
the components of n° = yJ\P\B°' satisfy 



(49) 



e e+ process 



Here, me is the electron mass: the 7 
will have a lower threshold than any other photon decay 
channel. Once again, high-energy photons from various 
directions can then constrain us to a finite region of n'^- 
space. 

The rate of energy loss for vacuum Cerenkov processes 
has been calculated by Altschul [18]. In particular, a 
charged particle with energy just above threshold [Ep — 
Eco + AE, with AE ^ Eco) will emit a photon with 
E~^ > AE with a decay rate of 



r = aZ^Try 



(AEy 

2EI 



(50) 



where a is the fine structure constant, Z is the particle's 
charge, and m is its rest mass. The mean free path of 
such a particle (assuming it to be moving with velocity 
V « Cp) can then be estimated as ^ = Cp/F. At higher 
energies, the charged particle will mainly lose energy to 
larger numbers of lower-energy photons, rather than a 
single photon that brings it below threshold. This pro- 
cess causes the energy to decrease even more rapidly; we 
should therefore view the estimate i above as an upper 
bound on the mean free path for a particle with energy 
Ep. 

In the case of photon decay, an exact expression for 
the lifetime of the photon is not yet known. However, 
we can estimate (see the Appendix) that the decay rate 
for photons above threshold energy E^o will be on the 
order of magnitude of 



F - aPB^E^ 



E^ 

E^ 

7 



(51) 



where a is the fine structure constant and i3 is a quantity 
of the same order of magnitude as the components of _B°. 
Using this estimate, one can then calculate a mean free 
path for photons as in the case of vacuum Cerenkov radi- 
ation. Roughly speaking, a photon well above threshold 
will have a mean free path of order [ajiB^)^^ times its 
Compton wavelength; if the photon is only barely above 
threshold, with energy E^ = (1 + e)Edo, its mean free 
path (in Compton wavelengths) is reduced by a factor 
of approximately y/e. 
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IV. EXPERIMENTAL CONSTRAINTS 

In general, the most stringent limits on the compo- 
nents of (kp)"''"''^ are those that arise from searches for 
vacuum birefringence [T]. However, the vanishing of 
the matrices (ke+Y-' 



(37a I implies that in 



and {ko-Y^ 

our case, bumblebee photons do not experience vacuum 
birefringence. We must thus turn to other experimental 
means of searching for Lorentz violation in the photon 
sector. In the following subsections, we will discuss lim- 
its arising from rotating electromagnetic resonator ex- 
periments, particle accelerator experiments, and cosmic- 
ray observations. 



A. Resonator experiments 

If the Maxwell field is Lorcntz-invariant, the frequen- 
cies of its modes in a resonant cavity will be independent 
of the cavity's orientation in space. However, if the pho- 
ton field propagates at different speeds in different di- 
rections, it is not hard to see that the resonant frequen- 
cies of the cavity can change if the cavity's orientation 
changes: this frequency depends on the "speed of light" 
in the cavity, and this speed is direction-dependent in 
our model if /? ^ and B 0. The magnitude of this 
frequency shift for a given set of matrices Ke± and ko± 
was calculated for a general cavity mode and geometry 
in |16j . In practise, this frequency shift is usually mea- 
sured by setting up two identical cavities, oriented at 
right angles to each other, and rotating these two cav- 
ities together. A difference in the speed of light within 
the plane of rotation would then show up as a "beat" 
between the frequencies of the two cavities, modulat- 
ing at twice the frequency of rotation. By looking for 
these "beats" at various points in the Earth's rotation 
and revolution, all eight independent components of the 
SME matrices {k^-Y'' and (ko+)'^ can in principle be 
measured. 

Sensitive measurements of the components of {ke-Y'' 
and {ko+Y'' have been performed by Herrmann et al. |19| 
and by Eisele, Nevsky and Schiller [5^. Both groups 
have bounded the components of {ke-Y'' to be C(10~^^) 
or less, and those of {ko+Y^ to be 0(10"-^^) or less. We 
can translate these measured bounds into a rough esti- 



mate of the bounds on our parameters (3 and i?^; given 



the dependencies given in ( 37b[ ) and ([37c 
expect the magnitude of nP = 



we would 
to be bounded 
below approximately 3 x 10"*, and the components of n 
to be approximately 3 x 10~^ or less.^ 

It is important to note, however, that both of the 



^ Note that the parameters of our model, as defined, are degen- 
erate; by rescaling our definition of B°' as above, we can set /3 
to ±1. Our physical parameter space is thus four-dimensional, 
with an additional discrete parameter (the sign of /9.) 





Best Fit 


la Conf. 


Intervals 




/3 > 


/3 < 


/3 > 


^ < 


n° 


0.69 


-1.61 


[-2.16,3.54] 


[-4.43,1.21] 




4.00 


-2.06 


[0.50,6.26] 


[-4.21,0.52] 




1.30 


3.61 


[-1.57,3.83] 


[0.07,5.76] 




-2.39 


-1.46 


[-5.29,0.52] 


[-3.88,0.96] 



TABLE I: Best fits and approximate la confidence intervals 
for (times 10~*) and n (times 10~*) in the standard Sun- 
centred frame, derived from the data of Stanwix et al. [21j . 



above mentioned groups [191 120j derived their respective 
bounds from their experimental data under the assump- 
tion that all eight independent components of (Ke_)*^ 
and (ko+Y'' could be varied independently. In our model, 
this is not the case; rather, as noted above, we have a 
four-dimensional parameter space (along with an addi- 
tional discrete parameter). A more thorough analysis 
should thus involve a (non-linear) regression on this pa- 
rameter space, starting from the experimental data. 

To perform such an analysis, we turn to the data of 
Stanwix et al. In Table I of [21 , the time- variation of the 
above-mentioned "beat" amplitudes are given. Applying 
a non-linear regression to these amplitudes gives the re- 
sults shown in Table|I]for the components in the Sun- 
centred frame. The best fit is found to occur for /3 > 0; 
however, a region of parameter space with /? < also 
falls in the overall la confidence region. The approxi- 
mate dimensions of the one-sigma confidence contour in 
both the /3 > and (3 < regions of parameter space are 
also given in Table |l] The point n'' = (0, 0, 0, 0) lies ap- 
proximately on the 75% confidence contour. While this 
might seem suggestive of a non-zero Lorentz-violating 
effect, such a confidence level can hardly be thought 
of as conclusive, especially considering that Stanwix et 
al. viewed as spurious the 2cr and 3cr signals found via 
their original analysis (with SME coefhcients assumed 
to be independent.) We must therefore conclude that 
resonator experiments have not yet observed a signal 
compatible with the generalized bumblebee model. 



B. Accelerator physics 

As noted above, vacuum Cerenkov radiation and pho- 
ton decay, though kinematically forbidden in conven- 
tional theories, are both allowed above a certain en- 
ergy threshold in the presence of Lorentz invariance in 
the photon sector. In the bumblebee photon case, this 
threshold will decrease as the components of the rescaled 
vector n" = -^/j/^fi?" increase. Direct observations of 
high-energy particles can therefore help us constrain our 
theory. 

An analysis along these lines has been performed by 
Hohensee et al. [22], through analysis of the operation 
of the LEP experiment and the Tevatron. In the case of 
vacuum Cerenkov radiation, they note that a threshold 
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energy Eco more than a few MeV below the electron 
and positron beam energies at LEP (i?LEP = 104.5 GeV) 
would have caused beam energy losses significant enough 
to be immediately apparent. In our case, this implies 
thafio 



k\ < 



4.9 X 10~ 



LEP 



(52) 



Moreover, from the known bounds due to resonator ex- 
periments (above), we know that the components of n 
are of order 10~^. This therefore implies that in the case 
/3 < 0, we must have |n°| < 4.9 x 10~^. 

For the case of photon decay, Hohensee et al. note that 
a significant fraction of predicted high-energy photons 
[E^ > 300 GeV) have been observed in the DO detector 
at Fermilab (specifically, in the study of isolated-photon 
production with an associated jet.) This implies that the 
threshold energy for photon decay cannot greatly exceed 
300 GeV; again using the bounds on the components of 
n from resonator experiments, we find that in the case 
/3 > we must have 



A < 



2me 



3.4 X 10" 



(53) 



More recently, work by Altschul [53] has extended 
these bounds by examining the amount of synchrotron 
radiation observed at LEP. The argument proceeds sim- 
ilarly to the analysis of vacuum Cerenkov radiation in 
LEP, given above. In a Lorentz-violating theory, the 
velocity of a relativistic charged particle moving in a 
uniform magnetic field will deviate from the "expected" 
velocity (i.e., that in the absence of Lorentz violation) 
by dv — k, where k is a linear combination (with coeffi- 
cients of order unity) of the Lorentz-violating coefficients 
{ke-Y^ , {ko+), and Ktr f24|. The fractional deviation of 
the synchrotron power radiated from its expected value 
is then given by of order 



AP 



87'^, 



(54) 



where 7 is the boost factor of the relativistic particle. 
When we integrate this power deviation over a full cycle 
of the charged particle, the parity-odd coefficients must 
drop out due to symmetry; thus, only {ke±y^ and Ktr 
can in principle contribute to the integrated synchrotron 
power loss. 

The greatest boost factor 7 of particles achieved at 
LEP was above 2 x 10^; the integrated power loss of 
particles in the storage rings at LEP was measured to 



The propagation direction k is, of course, not a constant for the 
electrons in a circular accelerator. However, the independent 
bounds on n from resonator experiments make this considera- 
tion moot. 



agree with the predictions of standard (non-Lorentz- 
violating) electrodynamics to within a fractional preci- 
sion of \AP\/P < 2 X 10"-* [ig. We can thus conclude 
that 



k < 



1 |AP| 
872 P 



6 X 10' 



-16 



(55) 



Since resonator experiments bound the components of 
(ke-)*"' to an order of magnitude below this, this bound 
is therefore only a bound on k^,-. In our case, ktr is given 
by equation ( 37d ). We can thus conclude that the known 



bounds on synchrotron radiation at LEP limit only the 
value of kti-] specifically. 



< ^/2\ktr\ < 3.5 X 10"^ 



(56) 



Note that this is a two-sided bound: since LEP is sen- 
sitive to AP/P being either positive or negative, and 
since the sign of AP/P is dependent on Ktr, then LEP 
measurements bound Ktr both above and below zero. In 
the current model, this means that nP is bounded both 
for ^ > and ^ < 0. 

This bound is comparable to the current bounds on 
the spatial components of n obtainable from resonator 
experiments. However, this bound should be taken cor- 
rect only to within an order of magnitude, due to our 
lack of knowledge about the precise functional form of 
k in our theory. A more precise estimate would involve 
a full calculation of the rate of synchrotron radiation in 
our model; however, such a calculation is outside the 
scope of this paper. 



C. Cosmic ray observations 

While man-made particle accelerators can impart very 
high energies to particles, it is well-known that "natural 
particle accelerators" elsewhere in the Universe put these 
efforts to shame; cosmic rays have been observed with 
energies of up to eight orders of magnitude more energy 
than the most energetic particles created in the labora- 
tory to date. Since our bounds on the parameters of our 
theory scale inversely with the energies of observed par- 
ticles, we will find that subject to some caveats (detailed 
below) cosmic ray observations give us the best bounds 
on the components of n'' . 

In the case of vacuum Cerenkov radiation, the Pierre 
Auger Collaboration has observed a few dozen cosmic 
ray showers with total energies above 57 EeV [23^. It 
is generally assumed that the primary particles for such 
cosmic ray showers are hyperrelativistic protons or nu- 



clei. In our case, the inequality (48) becomes more re- 



strictive as the mass of the primary particle decreases. 
Thus, if we wish to place conservative bounds on the 
parameters of the theory, we should assume that the 
primary is a fairly heavy nucleus. It is unlikely, how- 
ever, that these primaries are significantly heavier than 
^^Fe. We therefore assume that all the events listed in 
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Cosmic rays 


VHE photons 


|n"| 


6.8 X 10"^° 


2.1 X 10"* 


|n^| 


5.6 X 10"^" 


6.7 X 10"* 


\n^\ 


8.2 X 10"^" 


3.3 X 10"* 


W\ 


14.6 X 10"^" 


5.3 X 10"* 




19.9 X 10"^^ 
-2.3 X 10"^" 


46.6 X 10"^'' 
-3.3 X 10"^*^ 




6.0 X 10"^'^ 


2.7 X 10"^° 



TABLE II: Bounds placed on the components of = 
^/\j3\B°^ in the Sun-centred frame by cosmic-ray observations 
(/3 < 0) and high-energy gamma-ray observations (/3 > 0). 



[25] have a primary mass of MpCp — 52.1 GeV. If these 
events were to later be discovered to have a smaller mass 
M'p, it would simply rescale our limits on n'^ by a factor 
ofM^/Mp. 

The 27 events listed in [25_ thus disallow large vol- 
umes of parameter space. The remaining, allowed region 
of parameter space is a complicated non-uniform poly- 
choron, symmetric with respect to reflection about the 
origin. The maximum allowed magnitudes of each of the 
components of n'' are listed in Table |llj along with the 
maximum and minimum allowed values of the norm of 
n", and the volume Vn of the allowed polychoron. 

In the case (3 > 0, we would expect to see a (possi- 
bly direction-dependent) cutoff in the very high-energy 
(VHE) gamma-ray spectrum. Over the past decade, 
gamma-ray observatories such as HESS, VERITAS and 
MAGIC have catalogued several dozen VHE gamma- 
ray sources, and in many cases have been able to as- 
sociate these sources with known objects (either in the 
Milky Way or extragalactic.) However, two problems 
arise when attempting to use these sources to bound the 
parameters of our theory. The first is simply a mat- 
ter of orders of magnitude; the ratio of the energy of 
these gamma-rays to the electron mass is two orders of 
magnitude smaller than the boost factor of the charged 
particles that make up cosmic rays. The highest-energy 
photons detected thus far have energies of order 75 TeV, 
yielding a bounding factor 2mf,Cp/ Ej of approximately 
10"*; by contrast, the inverses of the boost factors for 
the charged-particle gamma rays observed by the Pierre 
Auger Observatory are of order 10"^*^. Thus, we cannot 
expect nearly so tight a bound on the components of 
as we obtained in the case /3 < 0. 

Second, the vast majority of the sources we can 
reliably use for this purpose lie in the plane of the 
Milky Way. Most known extragalactic sources of VHE 
gamma rays are associated with blazars located outside 
of the Local Supercluster, with non-negligible cosmolog- 
ical redshift (z > 0.05.) Once cosmological effects be- 
come non-negligible, our assumptions of nearly-flat met- 
ric ([s]) and nearly-constant vector field ^ can no longer 
be expected to hold. Instead, we would expect that 
would be approximate a Friedman-Robertson- Walker so- 



Name 


RA (°) 


5{ 


') Smax (TeV) ReL 


Crab Nebula 


83 


6 


22 


01 


75 


M 


RX J1713.7-3946 


258 


4 


-39 


76 


47 




Vela X 


128 


8 


-45 


60 


45 


M 


HESS J1825-137 


276 


5 


-13 


76 


40 


M 


MSH 15-52 


228 


5 


-59 


16 


40 


M 


Galactic Centre 


266 


3 


-29 


00 


32 


131] 


HESS J1809-193 


272 


6 


-19 


30 


30 


[32] 


HESS J1708-443 


257 





-44 


35 


20 


M 


LS 5039 


276 


6 


-14 


85 


20 


M 


RCW 86 


220 


7 


-62 


45 


20 


M 


HESS J1616-508 


244 


1 


-50 


90 


20 


M 


HESS J1813-178 


273 


4 


-17 


84 


20 


M 


M87 


187 


7 


12 


39 


20 


M 


Westerlund 2 


155 


8 


-57 


76 


18 


m 


HESS J1837-069 


279 


4 


-6 


95 


15 


M 


Kookaburra 


214 


5 


-60 


98 


15 


M 


HESS J1718-385 


259 


5 


-38 


55 


15 


M 


RX J0852.0-4622 


133 





-46 


37 


10 


m 


Cassiopeia A 


350 


9 


58 


82 


10 


M 


HESS J1702-420 


255 


7 


-42 


02 


10 


m 


HESS J1804-216 


271 


1 


-21 


70 


10 


m 


CTB 37A 


258 


6 


-38 


57 


10 


m 


Centaurus A 
AE Aquarii 


201.4 
310.0 


43.02 
-0.87 


5 

3.5 


m 
m 



TABLE HI: VHE gamma-ray sources with significant ob- 
served luminosity over 10 TeV. Shown are the names of the 
sources, their right ascension and declination (in degrees), 
and -Bmax, the approximate maximum energy of gamma rays 
observed in their spectra. The sources AE Aquarii and Cen- 
taurus A have -Emax < 10 TeV but are included due to their 
higher galactic latitude (—24.42° and 19.42°, respectively); 
see text. 



lution, and that the value of i?" would be subject to 
cosmological evolution. This would cause significant de- 



viations from the predictions of Section HI C which were 
predicated on ry"'' being constant throughout the propa- 
gation of the photon. For the purposes of this paper, we 
will therefore limit ourselves to VHE gamma-ray sources 
known or suspected to be within the Local Supercluster. 

With these caveats in mind, we proceed. A sur- 
vey of the literature reveals 22 known VHE gamma-ray 
sources with sig nific ant flux at energies greater than 10 
TeV (see Table [ml)" Only one of these sources, M87 
|36j . is extragalactic (with b = 74.49°). To improve 
our bounds in the directions orthogonal to the galac- 
tic plane, we therefore include two lower-energy sources: 



The large majority of the energies listed from these tables are 
extracted from plots in the references; they are therefore neces- 
sarily somewhat imprecise. 
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the radio galaxy Centaurus A [32] , with galactic latitude 
b = —24.42°, and the cataclysmic variable AE Aquarii 
|43j . with b = 19.42°. The resulting bounds are shown 
in Table [ill expected, they are significantly less strin- 
gent than the bounds imposed by charged-particle cos- 
mic ray observations; however, they are still competi- 
tive with the limits imposed on the components of ft by 
laboratory and accelerator experiments, and are more 
stringent in the case of n". 



V. DISCUSSION 

Laboratory experiments (both electromagnetic res- 
onators and bounds from accelerator physics) limit the 
components of to the 10~® level. In the case 
/3 < 0, cosmic-ray observations push these bounds down 
0(10^^") for all components of n^; for /3 > 0, similar ob- 
servations bound all the components of to the 10~® 
level. 

A few notes concerning the above bounds are in order. 
First, the bounds extracted from the regression in Sec- 



tion IV A are not based on the current strongest bounds 
on Lorentz violation in the photon sector, but rather on 
older work with more detailed reporting of data. One 
would expect that if the data underlying the most sen- 
sitive experiments to date [13 [3D] were analyzed in this 
way, our sensitivity to the components of n'^ would in- 
crease by approximately a factor of three (i.e., half an 
order of magnitude.) 

The great majority of the cosmic-ray sources and 
events used in Section IIV CI are in the Southern Celes- 
tial Hemisphere. This is simply due to the locations of 
the observatories in question. In the case of high-energy 
charged particles, the Pierre Auger Observatory (located 
in Argentina) does not yet have a Northern-Hemisphere 
counterpart; Pierre Auger North, a planned counterpart 
in Colorado, will not see "first light" for some years. 
In the case of high-energy gamma rays, this asymme- 
try is simply due to the fact that the HESS telescope 
(located in Namibia) has been in operation longer than 
the similar MAGIC and VERITAS telescopes (located in 
the Canary Islands and Arizona, respectively.) As these 
latter two telescopes and the air-shower arrays MILA- 
GRO (New Mexico), AS7 (Tibet), and ARGO (Tibet) 
report more high-energy sources (especially in the re- 
gion of the galactic plane lying in the Northern Celes- 
tial Hemisphere) , we can expect these bounds to become 
more symmetric. 

Finally, it is important to note that there are funda- 
mental limits on how well gamma-ray experiments can 
bound the parameters of our theory. Photons with suf- 



Note that a pair of (hypothetical) antipodal sources do not place 
the same bounds on the components of n'^; one bounds + 
n ■ k\, while the other bounds — n ■ k\. 



ficiently high energies can interact with background ra- 
diation (particularly the cosmic microwave background) 
and produce electron-positron pairs [33] • In particular, 
the "gamma-ray horizon" for 100 TeV gamma rays 
is approximately the size of the Local Supercluster, and 
for ^ 1 PeV gamma rays is approximately the size of the 
Milky Way. This pair-production attenuation therefore 
places fundamental bounds on the region of parameter 
space that can be bounded by gamma-ray observations: 
we do not expect to get extremely high-energy gamma 
rays from outside our own galaxy, but our bounds on 
the component of n orthogonal to the galactic plane 
will generically be insensitive to sources from within the 
Galaxy. 
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Appendix: Photon Lifetime 

In using particle physics to place bounds on Lorentz- 
violating theories, it is informative to examine the mean 
free paths of particles or, equivalently, their lifetimes. 
The rate of energy loss to vacuum Cerenkov radiation 
in a general Lorentz-violating theory has been previ- 
ously derived [TH]. However, the decay 7 — >■ e+e", as 
occurs in our case when /3 > 0, is less well-understood. 
The full quantum theory of Lorentz-violating photons is 
not yet known, and so a calculation of the exact pho- 
ton decay rate is impossible. However, it is well-known 
that when certain Lorentz-violating coefficients are suf- 
ficiently small, they can be "moved" between sectors 
via redefinitions of the metric (or, equivalently, redef- 
initions of the coordinates.) In particular, for a the- 
ory of Lorentz-violating photons without vacuum bire- 
fringence on a flat metric, it is possible to shift the 
Lorentz-violating coefficients entirely into the electron 
sector instead [15^ . The quantum field theory of Lorentz- 
violating electrons being much better understood j45) . 
we pursue this tactic. We will see that such a technique 
will only yield decay rates that are accurate to first or- 
der in j3; however, this calculation will still be valuable 
in estimating the mean free paths of Lorentz-violating 
photons. 

A similar calculation for "isotropic Lorentz-violating 
photons" has been performed by Hohensee et al. j22j : 
we will follow their techniques here. For consistency 
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with this paper and f?5] , we use a metric with signature 
(+, — , — , — ). In particular, when switching to this sign 
convention, our definition of fj"''' will change; in this sign 
convention, a metric of the form 



f^ab 



V 



ah 



(A.l) 



is physically equivalent to the metric used in the previous 
sections. 

We start with the "Lorentz-violating QED" La- 
grangian for bumblebee photons minimally coupled to 
an electron field ^ in flat spacetime: 



1 



(A.2) 



with fj defined as in (A.l I. The process used by Ho- 



hensee et al. to derive the photon lifetime consists of 
several steps: 

1. We redefine the metric so that the "photon metric" 
the "true metric" of the theory. Since the 
gamma matrices in the electron kinetic term are 
defined with respect to the original metric -q"-^ (i.e., 
{t'^iT''} — 2?7°^), we must rewrite these matrices 
in terms of new gamma matrices 7° defined such 
that {7^,7''} = 27/°''. To first order in j3, these 
matrices are related by 



7 



(A.3) 



Applying this definition to our Lagrangian, we find 
that our new Lagrangian is 



FabF'''' 

.(A.4) 

where indices are now raised and lowered with the 
metric fj'^^. We will hereafter "drop the tildes" for 
notational convenience. 

2. The kinetic term for the electron now contains 
non-standard time derivatives; these must be elim- 
inated to successfully quantize the theory fTE\ . To 
do this, we define a new spinor field x such that 
"0 = Ax, where the matrix A satisfies 



7° + ?S°bS) A 



(A.5) 



(A.6) 



'7" I -7" + -ti"B""(h I A = 1. 
To first order in /3, this implies that 
A = l 

Rewriting the Lagrangian in terms of x then yields 

C = Ir'xiaDbX ~ mxx - lFabF-\ (A.7) 
where 



fjab ^^ab_l^ (^{B" f7f'' + 2B°B^''tj^^° - B" B''^ (A.. 



and 



/3 



(A.9) 



3. Write down the invariant matrix element Ai. Us- 
ing the photon and electron polarizations defined 
in |22j . we find this to be given by 



iMrs = -ieea{p)r''u'-'-HqhbV^'\k) 



(A.IO) 



where p"" is the incoming photon momentum, (p) 
is its polarization, g° and k"' are the outgoing 
electron and positron momenta respectively, and 
u^''^ (q) and w^*^ (k) are their respective polarization 
states. To obtain the photon lifetime, we will need 
to sum over the final fermion polarization states 
and average over photon polarization states; using 
the trace identities derived in [IS], this yields 



r,s I 
e r.s 

= 2e^{r''fiab{m^ + fi%fi,aq'k'') 

~2fia%br''v''qdke). (A.ll) 

We have used the spinor normalization conventions 
chosen in [22] (i.e., N{q) = 2Eq and similarly for 
fc.) 

4. Integrate over the final electron and positron mo- 
menta to obtain the photon lifetime as a function 
of its energy E-y and propagation direction p: 



m 



1 1 



d'V d^fc 
2Eq 2Ek 



(A.12) 



I have as yet been unable to derive a closed-form an- 
alytical expression for this lifetime. However, an order- 
of-magnitude estimate of the photon decay rate may be 
obtained by estimating the quantity \M.\'^ / AEqE^ for an 
on-shell decay and multiplying this quantity by the al- 
lowed volume of phase space; in other words, 

j3; 



dV d^fc 
2Eq 2Ek 



IMl^S^pf" - qf" - kf") 



X / d^qd^kS'^{p^' - - k^') 



4.EqEk 



This equation can be thought of as defining the quantity 
in angle brackets above. To estimate its order of magni- 
tude, we can use the dispersion relation for the electrons 
in this theory. 



fiac'nb''k°'k^ = 



(A.13) 
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along with momentum conservation, p° = k"' + q"", to 
put (A.lll in the form 



|X|2 = 2e2 ^2(3E'^ (^-2B"B-p+{B-pf' 

+ 2(1 - /3(5")2)m2 + /3(B,r )2 + /3{Baq^f 

+ (A.14) 

All of these terms except for the last two are constant 
over the mass shell. To estimate the order of magni- 
tude of these last two terms, we note that the electron 
dispersion relation satisfies i?^ — P = + 0{P)] thus, 



a\2 



fi{B°Ek - |B||fc|cos( 



\B\^|El 



(A.15) 



the equation above is of equal or lesser magnitude than 



the first, since pB^E'^i > pB^E, 
Thus, we can estimate that 



•) 



(A.18) 



where a is the fine structure constant. The volume of 
kinematically accessible phase space, meanwhile, can be 
shown to be 




m 



(A.19) 



where we have discarded the cos 9 term because it is of 
order unity. We can further estimate that for a generic 
decay E^ (and Eg) will be of order E^/2, and that (since 
E^ > Edo ^ nr^) the resulting electrons will have rela- 
tivistic velocities. This then implies that 



(A.16) 



This integral is most easily done by changing coordinates 
on phase space to u = q + k and v = ^{q — k); the re- 
sulting expression sets u = p, leaving an expression pro- 
portional to the volume of an ellipsoidal shell in w-space. 
We can then perform a further linear transformation on 
V to make this shell spherical. Combining these results, 
we can then estimate that the photon lifetime is 



where ;B is a quantity of the same order as the compo- 
nents of We can then estimate the value of |A^p to 
be 



-(3B^E^. 1 - 



(A.20) 



2e^ [PB^E?^ + 2(1 - l3{B° f)m^] 



26^/3^2^2 



7' 

where we have redefined B to include the contributions 



of the first term in (A.14 1. (Note that the second term in 



to within an order of magnitude. Note that the scaling 
behaviour of this result is consistent with the exact re- 
sults derived in the spatially isotropic case by Hohensee 
et al. [221 • 
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